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£\) ■ Shell-model states involving several pseudospin doublets and "intruder" levels in 

f^*^ ■ nuclei, are combined into larger multiplets. The corresponding single-particle spec- 

' trum exhibits a supersymmetric pattern whose origin can be traced to the rela- 

f**"^ | tivistic pseudospin symmetry of a nuclear mean-field Dirac Hamiltonian with scalar 

■^J- , and vector potentials. 

o 



1. Introduction 



O ' Pseudospin doublets 1 in nuclei refer to the empirical observation of quasi- 

degenerate pairs of certain shell-model orbitals with non-relativistic single- 
nucleon radial, orbital, and total angular momentum quantum numbers: 



(n,£,j =£+1/2) and (n - l,t + 2,j' = t + 3/2) . (1) 



^ | The doublet structure (for n > 1) is expressed in terms of a "pseudo" 

orbital angular momentum £ — t + 1 and "pseudo" spin, s = 1/2, coupled 
to j = £ — 1/2 and j' = £ + 1/2. For example, [nSi/ 2 , (n — 1)g?3/2] will 
have I = 1, etc. The states in Eq. (1) involve only normal-parity shell- 
model orbitals. The states (n = 0,1, j = I + 1/2), with aligned spin and 
no nodes, are not part of a doublet. This is empirically evident in heavy 
nuclei, where such states with large j, i.e., Ogg/2, 0/in/2, 0ii3/2, are the 
"intruder" abnormal-parity states, which are unique in the major shell. 

Pseudospin symmetry is experimentally well corroborated and plays a 
central role in explaining features of nuclei 2 including superdeformation 3 
and identical bands 4 . It has been recently shown to result from a relativistic 
symmetry of the Dirac Hamiltonian in which the sum of the scalar and 
vector nuclear mean field potentials cancel 5 . The symmetry generators 6 
combined with known properties of Dirac bound states, provide a natural 
explanation 7 for the structure of pseudospin doublets and for the special 
status of "intruder" levels in nuclei. 



1 
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Figure 1. Nuclear single-particle spectrum Figure 2. Typical supersymmetric pat- 
composed of pseudospin doublets and an tern. The Hamiltonians Hi and H2 have 
"intruder" level. All states share a common identical spectra with an additional level 
i and s = 1/2. The corresponding Dirac k- for Hi when SUSY is exact. The opera- 
quantum numbers are also indicated. tors L and L' connect the partner states. 



Figure 1 portrays the level scheme of an ensemble of pseudospin dou- 
blets, Eq. (1), with fixed £, j, j' and n = 1,2,3,... together with the 
"intruder" level (n = 0, £,j = I + 1/2). The single-particle spectrum ex- 
hibits towers of pair-wise degenerate states, sharing a common £, and an 
additional non-degenerate nodeless "intruder" state at the bottom of the 
spin-aligned tower. A comparison with Fig. 2 reveals a striking similarity 
with a supersymmetric pattern. In the present contribution we identify 8 
the underlying supersymmetric structure associated with a Dirac Hamilto- 
nian possessing a relativistic pseudospin symmetry. 



2. Supersymmetric Quantum Mechanics 

Supersymmetric quantum mechanics (SUSYQM), initially proposed as a 
model for supersymmetry (SUSY) breaking in field theory 9 , has by now 
developed into a field in its own right, with applications in diverse areas of 
physics 10 . The essential ingredients of SUSYQM are the supersymmetric 
charges and Hamiltonian 

which generate the supersymmetric algebra 

[H,Q±] = {Q±,Q±} = , {Q_,Q+} = W. (3) 

The partner Hamiltonians Hi and H2 satisfy an intertwining relation, 
LH\ = H2L, where in one-dimension the transformation operator L = 
+ W{x) is a first-order Darboux transformation expressed in terms of 
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a superpotential W(x). The intertwining relation ensures that if is an 
eigenstate of Hi, then also $ 2 = L^i is an eigenstate of H 2 with the same 
energy, unless L^i vanishes or produces an unphysical state (e.g. non- 
normalizable). Consequently, as shown in Fig. 2, the SUSY partner Hamil- 
tonians Hi and H2 are isospectral in the sense that their spectra consist 
of pair-wise degenerate levels with a possible non-degenerate single state in 
one sector (when the supersymmetry is exact). The wave functions of the 
degenerate levels are simply related in terms of L and L^. Such character- 
istic features define a supcrsymmetric pattern. In what follows we show 8 
that a Dirac Hamiltonian with pscudospin symmetry obeys an intertwining 
relation and consequently gives rise to a supersymmetric pattern. 

3. Dirac Hamiltonian with Central Fields 

A relativistic mean field description of nuclei employs a Dirac Hamiltonian, 
H = a-p + fi(M + Vs) + Vv, for a nucleon of mass M moving in external 
scalar, Vs, and vector, VV, potentials. When the potentials are spherically 
symmetric: Vs = Vs(r), Vy = Vv(r), the operator K = —(3(cr • I + 1), 
(with cr the Pauli matrices and £ = —ir x V), commutes with H and its 
non-zero integer eigenvalues k = ±(j + 1/2) are used to label the Dirac 
wave functions 



Here G K (r) and F K (r) are the radial wave functions of the upper and lower 
components respectively, Yi and \ are the spherical harmonic and spin 
function which are coupled to angular momentum j with projection m. 
The labels k = -(j + 1/2) < and £' = I + 1 hold for aligned spin 
j = £+1/2 (si/ 2 ,p 3 /2, etc.), while k = (j + 1/2) > and £' = I - 1 hold 
for unaligned spin j =1 — 1/2 (pi/2, d 3 / 2 , etc.). Denoting the pair of radial 
wave functions by 





(5) 



the radial Dirac equations can be cast in Hamiltonian form, 

H K $ K = E $ K , 



(6) 



with 




(7a) 



A(r) 



Vs + V v , £(r) = V s -W. 



(7b) 
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The nuclear single-particle spectrum is obtained from the valence bound- 
state solutions of Eq. (6) with positive binding energy (M—E) > and total 
energy E > 0. The non-relativistic shell-model wave functions are identified 
with the upper components of the Dirac wave functions (4) . For relativistic 
mean fields relevant to nuclei, V$ is attractive and Vy is repulsive with 

typical values V s (0) 400, V v (0) ~ 350, MeV. The potentials satisfy 

rVs, rVy — ► for r — > 0, and Vs, Vv — ► for r — > oo. Under such 
circumstances one can prove 7 the following properties which are relevant 
for understanding the nodal structure of pseudospin doublets and intruder 
levels in nuclei. 

(a) The radial nodes of F K {rip) and G K (no) are related. Specifically, 

n F = n G for k < , , . 

n-F = n G + 1 for k > . 

(b) Bound states with n F = tiq = can occur only for n < 0. 



4. Relativistic Pseudospin Symmetry in Nuclei 

A relativistic pseudospin symmetry occurs when the sum of the scalar and 
vector potentials is a constant 

A(r) = Vs(r)+lV(r)=Ao. (9) 

A Dirac Hamiltonian satisfying (9) has an invariant SU (2) algebra gener- 
ated by 6 

MU)- ,10) 

Here = cr^/2 are the usual spin operators, s M = U p s^U p and U p = 
is a momentum- helicity unitary operator 11 . In the symmetry limit the 
Dirac eigenfunctions belong to the spinor representation of SU(2). The 
relativistic pseudospin symmetry determines the form of the eigenfunctions 
in the doublet to be 



y K1 <o, m = z [ _ r „ L K! = -e<o, j = e-i/2, (iia) 



l(G Kl [Yi-iX]P 
1 / G K2 [Y i+1 xf [ 



* K2 >0,m = -| .„ r ;t l( J'1 I K 2-i?+l>0, j' =1+1/2, (lib) 



r 



iF K2 [YiX] 



and imposes the following conditions on their radial amplitudes: 
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Figure 3. Top left panel: the up- 
per components g(r) = rG K (r) of 
the 2S]/2, (solid line) and 1^3/2 
(dashed line) Dirac eigenfunctions 
in 208 Pb. Top right panel: testing 
the differential relation of Eq. (12b) 
for the upper components of 2s]y 2 
(ki = —1) and lrf 3/ / 2 («2 = 2). 
Bottom panel: the lower compo- 
nents f(r) = rF K (r) of 2sj/2 and 
l<^3/2> testing relation (12a). Based 
on calculations in [12,13]. 



r (Fermi) 



F K 



F K 



dG Kl ki„ _ dG K2 K2 n 

1 <-*Kl — 1 ^2 

ar r ar r 



(12a) 
(12b) 



The two eigenstates in the doublet are connected by the pseudospin gener- 
ators Sfj, (10). The lower components are connected by the usual spin oper- 
ators and, therefore, have the same spatial wave functions. Consequently, 
the two states of the doublet share a common I which is the orbital angu- 
lar momentum of the lower component. The Dirac structure then ensures 
that the orbital angular momentum of the upper components in Eq. (11) is 
I = 1-1 for j = 1-1/2 = £+1/2, and £+2 = £+1 for j' = £+1/2 = £+3/2. 
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This explains the particular angular momenta defining the pseudospin dou- 
blets in Eq. (1). The radial amplitudes of the lower components are equal 
(12a) and, in particular, have the same number of nodes np = n. Prop- 
erty (a) of the previous section then ensures that G Kl has n nodes and G K2 
has n — 1 nodes, in agreement with Eq. (1). Property (b) ensures that the 
Dirac state with rif = ng = 0, corresponding to the "intruder" shell- model 
state, has a wave function as in Eq. (11a) with n < 0, and does not have a 
partner eigenstate (with k > 0). 

Realistic mean fields in nuclei approximately satisfy condition (9) with 
Ao ~ 0. The required breaking of pseudospin symmetry in nuclei is small. 
Quasi-degenerate doublets of normal-parity states and abnormal-parity lev- 
els without a partner eigenstate persist in the spectra. The relations (12) 
between wave functions have been tested in numerous realistic mean field 
calculations in a variety of nuclei, and were found to be obeyed to a good 
approximation, especially for doublets near the Fermi surface 12 > 13 . A rep- 
resentative example for neutrons in 208 Pb is shown in Fig. 3. 

5. Relativistic Pseudospin Symmetry and SUSY 

In the pseudospin limit, Eq. (9), the two Dirac states W Kl <o,m an d ^ K2 >o,m 
of Eq. (11) with K\ + k 2 = 1 are degenerate, unless both the upper and 
lower components have no nodes, in which case only \& Kl <o,m is a bound 
state. Altogether, as shown in Fig. 4, the ensemble of Dirac states with 
Ki + «2 = 1 exhibits a supcrsymmetric pattern of twin towers with pair- wise 
degenerate pseudospin doublets sharing a common £, and an additional non- 
degenerate nodeless state at the bottom of the K\ < tower. An exception 
to this rule is the tower with k 2 = 1 {pi/2 states with I = 0), which is 
on its own, because states with K\ = do not exist. The supersymmetric 
structure arises because in the pseudospin limit the respective radial Dirac 
Hamiltonians, H Kl and H K2 , satisfy an intertwining relation of the form 



LH Kl = H K2 L , 



(13) 



with «i + «2 = 1. The transformation operator is found to be 




(14) 



L connects the two doublet states 



L$ K1 =b(M + A -E)$, 



(«i + k 2 = 1) , 



(15) 
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Figure 4. Schematic supcrsymmetric pattern in the pseudospin limit of the Dirac Hamil- 
tonian [8]. 

and identifies the two states as supersymmetric partners. Eq. (15) relies 
on the input that the two states are eigenstates of the Dirac Hamiltonian 
with E Kl = E K2 = E and their wave functions satisfy the relations in 
Eq. (12). The fact that nuclear wave functions, obtained in realistic mean- 
field calculations, obey these relations to a good approximation, confirms 
the relevance of supersymmetry to these nuclear states. 

Constructing supersymmetric charges Q± and Hamiltonian TL from L 
and H Kl , H K2 as in Eq. (2), ensures the fulfillment of Eq. (3), except for 
the last relation which now reads 

{Q-,Q+} = b 2 [H - (M + A )][H - (M + A )] . (16) 

Eq. (16) involves a polynomial of TL, indicating a quadratic deformation 
of the conventional supersymmetric algebra 14 . The latter arises because 
both the Dirac Hamiltonian, iJ re , and the transformation operator, L, are 
of first order. In real nuclei, the relativistic pseudospin symmetry is slightly 
broken, implying A(r) ^ A in Eq. (9). Taking H K as in Eq. (7) and L as 
in Eq. (14) but with A — ► A(r), we now find that 

LH Kl -H K2 L = ib^-a 2 (17) 

Furthermore, {Q-,Q+} has the same formal form as in Eq. (16), but the 
appearance of A(r) instead of Ao implies that the anticommutator is no 
longer just a polynomial of TL. 
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6. Summary 

We discussed a possible grouping of shell-model single-particle states into 
larger multiplcts, exhibiting a supersymmetric pattern. The multiplets in- 
volve several quasi-degenerate pseudospin doublets and intruder levels with- 
out a partner eigenstate. In contrast to previous studies of pseudospin in 
nuclei, the suggested grouping of nuclear states treats the intruder levels 
and pseudospin doublets on equal footing. The underlying supersymmetric 
structure is linked with an approximate relativistic pseudospin symmetry 
of the nuclear mean-field Dirac Hamiltonian. The relativistic pseudospin 
symmetry imposes relations between the upper and lower components of 
the two Dirac states forming the doublet. These relations, which are obeyed 
to a good approximation by realistic mean-field wave functions, imply that 
the Dirac Hamiltonian with pseudospin symmetry obeys an intertwining 
relation which gives rise to the indicated supersymmetric pattern. 

Acknowledgments 

This work was partly supported by the U.S. -Israel Binational Science Foun- 
dation, in collaboration with J.N. Ginocchio (LANL), and partly by the 
Israel Science Foundation. 

References 

1. K.T. Hecht and A. Adler, Nucl. Phys. A137, 129 (1969); A. Arima, M. Harvey 
and K. Shimizu, Phys. Lett. B30, 517 (1969). 

2. A. Bohr, I. Hamamoto and B.R. Mottelson, Phys. Scr. 26, 267 (1982). 

3. J. Dudek et al, Phys. Rev. Lett. 59, 1405 (1987). 

4. W. Nazarewicz et al, Phys. Rev. Lett. 64, 1654 (1990); F.S. Stephens at al., 
Phys. Rev. Lett. 65, 301 (1990). 

5. J.N. Ginocchio, Phys. Rev. Lett. 78, 436 (1997). 

6. J.N. Ginocchio and A. Leviatan, Phys. Lett. B425, 1 (1998). 

7. A. Leviatan and J.N. Ginocchio, Phys. Lett. B518, 214 (2001). 

8. A, Leviatan, Phys. Rev. Lett. 92, 202501 (2004). 

9. E. Witten, Nucl. Phys. B188, 513 (1981). 

10. G. Junker, Supersymmetric Methods in Quantum and Statistical Physics, 
(Springer- Verlag, Berlin, Heidelberg 1996). 

11. A. L. Blokhin, C. Bahri, and J. P. Draayer, Phys. Rev. Lett. 74, 4149 (1995). 

12. J.N. Ginocchio and D.G Madland, Phys. Rev. C 57, 1167 (1998); GA. 
Lalazissis et al, ibid. 58, R45 (1998); J. Meng et al, ibid. 59, 154 (1999). 

13. J.N. Ginocchio and A. Leviatan, Phys. Rev. Lett. 87, 072502 (2001); J.N. 
Ginocchio, Phys. Rev. C 66, 064312 (2002). 

14. N. Debergh et al, J. Phys. A 35, 3279 (2002); L.M. Nieto et al, Ann. Phys. 
305, 151 (2003). 



